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ABSTRACT 

We propose a two-temperature radial inflow-outflow model near Sgr A* with self-consistent feeding 
and conduction. Stellar winds from individual stars are considered to find the rates of mass injection 
and energy injection. These source terms help to partially eliminate the boundary conditions on the 
inflow. Electron thermal conduction is crucial for inhibiting the accretion. Energy diffuses out from 
several gravitational radii, unbinding more gas at several arcseconds and limiting the accretion rate to 
< 1% of Bondi rate. We successfully fit the X-Ray surface brightness profile found from the extensive 
Chandra observations and reveal the X-Ray point source in the center. The super-resolution technique 
allows us to infer the presence and estimate the unabsorbed luminosity L rj 4- I0 32 erg s _1 of the point 
source. The employed relativistic heat capacity and direct heating of electrons naturally lead to low 
electron temperature T e ss 4 • I0 10 K near the black hole. Within the same model we fit 86 GHz 
optically thick emission and obtain the order of magnitude agreement of Faraday rotation measure, 
thus achieving a single accretion model suitable at all radii. 

Subject headings: accretion, accretion disks — conduction — Galaxy: center — stars: winds, outflows 



1. INTRODUCTION 



Our Galaxy hosts a supe rmassive black hole (BH) with 



a ma ss M = 4.5 • 10 b Af Q (Ghez et al. 



2008 



Reid et al. 



2008| at a distance R = 8A kpc. The Uti exhibits low 
luminosity state probably due to inefficient feeding and 
cooling. Almost all available matter outflows fr om the re- 
gion, whereas only the small fraction accretes ( Quataert 
2004|. This feeding region within several arcseconds 



contains X-Ray emitting gas, but some X-Rays are ex- 
pected from a synchrotron self-Compton (SSC) or syn- 
chrotron source from accretion at several Schwarzschild 
radii r s . The study of X-Rays offers a unique opportu- 
nity to test the full range of accretion scales from several 
" to r g = I0~ 5 " and construct a single model. 

Modeling the accretion flow with such a huge range of 
scales is a challenge. 3D SPH simulations are performed 
in the outer region between 1" and 10" ( Rockefeller et 
[all|2004[ [Cuadra et al.||2008| ). Latest MHD simulations 
( jsharma et al.|2U08 1 are limited to 3 orders of magnitude 
in radius an d axial symmetr y. Only the one-dimensional 
calculation (Quataert 2004) can in principle resolve the 
flow everywhere. Thus, ID modeling is the approach we 
adopt extending it down to the BH horizon. 



We analyze the quiescent observations (Muno et al. 
2008 1 of X-Ray emission from central several arcseconds 
around Sgr A* in ^J2j The total exposure is 25 times 
longer comp ared to previously analyzed data (Baganoff 



et al.||2003[ ). The super-resolution processing based on 
spacecraft dithering helps resolving sub-pixel scales. The 
up-to-date data on stellar wind emitters are summarized 
in Sj3] We smooth matter ejection rates of individual stars 
over radius and sum them into a single feeding rate, also 
properly averaging the wind velocity. This presents a sig- 
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nifican t improvement over an ad-hoc feeding in Quataert 
( 2004 ) . The dynamical two-temperature equations are 
derived in ^4j We consider the electron conduction the 
main energy transport mechanism, approximating the 
unsat urated heat flux by a simple formula. The Bondi 
flow ( |Bondi|[l952 ) without heat transport overestimates 
the X-Ray luminosity by a factor of I0 3 . The other im- 
portant effects considered are the relativistic heat capac- 
ity of electrons and superadiabatic heating equivalent to 
entropy production. The ways to solve the resulting sys- 
tem of equations and corresponding results are presented 
in S|5] We employ the shooting method and find the min- 
imum \ 2 fit for X-Ray surface brightness profile, simulta- 
neously fitting 86 GHz flux. The best fit model requires 
X-Ray point source. The viability of a non-cooling radial 
flow is examined. 

2. OBSERVATIONS 

Central several arcseconds of the Galaxy were observed 
quite often over the past several years. The rich re- 
gion contains poin t sources identified as X-Ray b inaries 
(|Muno et al. 2009 ) and extended emission features ( |Muno| 
et al.|2008[ ) together with the source coincident with Sgr 
A*. The latter is expected from hot acc reting gas, and 



sourc e confusion is practically impossible ( jBaganoff et al 
2003 1 . Sgr A* source exhibits signif icant X-Ray flares as" 



sociated with the SSC mechanism (Baganoff et al.|2001 | 
or synchrotron ( |Dodds-Eden et aT 2009) T We are in- 
terested in quiescent emission, so we exclude the flaring 
state. We bin the observations in 628 seconds as a com- 
promise between the time resolution and the number of 
counts. About 4 photons on average are received dur- 
ing 628 seconds and we take only the observations with 
less than 15 photons, thereby accumulating 953 ks in 
the quiescent state. The quiescent state also produces 
some point source X-Rays , likely associated with SSC 
HMoscibrodzka et al.||2009| . We model these by a PSF- 
broadened central point source. We eliminate the emis- 
sion from the point sources and bright extended sources 
offset from Sgr A* (see Figure [I]). The bright extended 
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Fig. 1. — Chandra image of central 6" around Sgr A*, 
sources and strong extended features are subtracted. 



Point 



emission may arise from the colliding winds of two strong 
close emitters or from the collision of hot outflowing ma- 
terial with cold molecular material. We exclude both 
effects from modeling of an averaged flow pattern. 

We construct the surface brightness profile in counts 
per pixel squared for the duration of observation as a 
function of distance from the BH. The size of Chandra 
pixel is 0.5", which may seem to pose a limit on radial 
binning of brightness profile. However, the position of 
satellite is not steady over the duration of observations, 
but is findable with the 0.1" accuracy by comparing with 
the known positions of bright point sources. Then we can 
achieve 0.1" super-resolution accuracy in surface bright- 
ness profile from knowing the orientation of the detector 
pixels at any given time. The final profile is shown on 
Figure^ (error bars) together with the point-spread func- 
tion (FSF) (dashe d) found from the nearby point source 
J174540.9-290014 ( |Muno et al.|2009| ). The PSF is scaled 
to match the contribution from the point source. The 
counts cease to be monotonic at about 5" due probably 
to the production of X-Rays in collisions of cold and hot 
regions. Therefore, only radiation within the central 5" 
is to be modeled. As we are interested in how symmetric 
the surface brightness profile is, we divide the emitting 
region into 4 sectors 90 deg each centered on Sgr A* and 
extract the surface brightness profile in each sector. The 
standard deviation of counts between sectors is below 2a 
the noise within 5", but rises to several a outward from 
5". This justifies our choice of the outer radiation bound- 
ary and proves the applicability of the radial model. Let 
us now look in more details on manufacturing of the X- 
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Fig. 2. — Observed radial surface brightness profile (error bars), 
best fit (solid) and the point source contribution to emission 
(dashed). The point source contribution is the scaled PSF. 



3. STELLAR WINDS FEEDING 

The Galactic Center region has a concentration of mas- 
sive Wolf-Rayet and blue giant stars, expellin g strong 
winds from their surfaces ( Martins et al.| 2007). As the 
strongest wind emitters are usually the brightest stars, 
all wind emitters are easily identifiable. We take the lat 



est data on ejection rate s and velocities (Martins et al. 
2007 Cuadra et al.|2008 ) an d complement them w i th the 



or bital p arameters o t stars (Paumard et al. 2006[ Lu et 
all[2009|. Following [Cuadra et al.| |2008f 



we minimize 

eccentricities for the stars not belon ging to the stellar 
disks as identified by Lu et al. ( 2009| . The w ind speeds 
v w and ej ection rates are taken directly from Cuadra et 
[all pOOl. 



er e are several ways to treat the winds. Rockefeller 
et al. ( 2004 1 perform ed a simulation with winds from 
steady stars, whereas Cuadra et al. (2008) considered 
moving stars. In both cases the time to reach the quasi- 
steady solution 300 — 1000 yrs is comparable to or longer 
than the orbital period at the stagnation point 350 yrs. 
Thus, it is reasonable to average over stellar orbits in 
a search for a steady-state prescription of feeding. We 
reconstruct the full 3D orbits, but retain only the apoc- 
enter and pericenter distances for the stars. We smooth 
the total wind ejection rate for each star over the ra- 
dial extent of its orbit and then smooth with the narrow 
Gaussian filter to eliminate the divergences at the turn- 
ing points. We add the resultant feeding profiles together 
to obtain the total feeding rate as a function of radius (see 
Figure [3| . We square average the wind velocities weigh- 
ing the contribution of each star by its mass loss rate. 
However, the winds also acquire the velocity of a star as 
viewed by a distant observer. We neglect stars' proper 
motions in calculations of wind energy. They are negli- 
gible at several arcseconds, but would rather contribute 
to the angular velocity of matter within 1", where feed- 
ing is dominated by few stars. The dependence of the 
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Fig. 3. — Mass input into the feeding region around the BH on the 
upper panel. Square averaged wind velocity v w on the lower panel. 
Feeding is averaged over stellar orbits. Each wiggle represents a 
turning point of a single orbit. Only S02 star feeds matter within 
0.8". 

averaged wind speed on radius is shown on Figure [3j 
Quataert (2004) assumed the power-law mass injection 
rate q(r) cx r~ 1 for r £ [2", 10"]. The power-law index 
r] = 2 corresponds to zero slope of M (r) cx r 2 q(r) (see 
Figure[3]) and agrees better with the present calculations, 
whereas their choice of constant wind velocity does not 
agree with the present estimate. 



We also incorporate S02 star (Martins et al. 20081 
into the calculations. The mass loss rate Mg 02 = 6 • 
lO _8 M0year _1 of S02 is taken to coincide with that of 
r Sco. S02 has a spectral type B — 2.51^ and a mass 



M « 16M ( |Mokiem et al ||2005[ |Martins et al.||2008| ), 
whereas r Sc o has a very close ty pe B0.2V and a mass 
M k, 15M ( |Mokiem et al. 1120051). The inferred accre- 
tion rate onto the black hole ( Sharma et al.|[2007a|b| 
3- 10 -8 MQyear -1 is actually smaller than Mso2, thus the 
whole accreted material can in principle be provided by 
a sing le weak wind emitte r. This result is very different 



from Cuadra et al. (2008), who assumed all the matter 
accretes from the inner boundary of the simulation, thus 
obtaining in a simplified treatment a much larger accre- 
tion rate. However, the direct feeding mechanism (Loeb 



2004 1 by S02 does not work, as its revised A/502 is much 
below the value required for feeding without the angu- 
lar momentum. In turn, th e direct feeding by IRS 13E3 
( Moscibrodzka et al.||2006| produces too large accretion 
rate in the absence or conduction. 

4. DYNAMICAL EQUATIONS 



4.1. Energy transport mechanism 

Radiatively inefficient flows can be mediated signif- 
icantly by th e energy transfer from the inner regions 
to the outer ( Blandford & Bcgclman 1999 Johnson & 
Quataert 2007; Sharma ct al. 2008]) . Such transfer hap- 
pens m two distinct ways: via convection or via diffusive 
energy transport. Convection is seen in n umerical simu- 
lation s. It happens via Alfven instability ( Igumenshchev 



2006[) and magneto-thermal instability (iVlllj ( |sharma 
et al.|2008j) and modifies the density profile. Let us show 



that the electron heat conduction wins over convection 
in the accretion flow. First, the MTI is driven by ther- 
mal conduction, at any moment the electron conduction 
flux is larger then the MTI-induced heat flux. Convec- 
tion implies the motion of large-scale magnetized eddies, 
which in turn split into smaller eddies and develop the 
whole turbulent cascade. In such settings t he electron 
conduction is on ly inhibited a factor of ~ 5 ( Narayan & 
Medvedev||2001|). The speed of electrons is a factor of 



■\Zm p /m e larger than the sound speed and the convec- 
tion is subsonic; the same factor lowers the ion diffusive 
heat transport. The relative strength of convective heat 
flux is proportional to the gradient of logarithmic en- 
tropy, which is normally weaker than the proportionality 
to the gradient of logarithmic temperature of conductive 
flux. Combining both effects wc conclude that, if there 
is convection or diffusion, then there is stronger conduc- 
tion. Severe inhibition of electron conduction happens, 
if the turbulent cascade does not develop and mixing is 
absent. This is not the case when the gas accretes. The 
strength of turbulent magnetic field increases then in the 
con vergent flow leading to dissipation and effective mix- 
ing ( Shvartsman||197l Shcherbakov||2008' ) . It is reason- 
able to - EHrlFtEiatThe wli^le - turT5uTent cascade develops 
and the electrons relatively freely find their way around 
magnetic field lines to connect the different regions of 
the flow. When the electrons and ions get decoupled 
from each other, the ion entropy may get equilibrated by 
convection, whereas the electron temperature levels due 
to conduction. The investigation of this possibility is left 
for future research. In present paper we take the energy 
transport to happen solely via electron conduction. 

There are several different regimes of conduction. 
First, the collisionality of the flow changes from the large 
radii to the inner radii as the mean free path of par- 
ticles I exceeds the flow size r. As the flow gets only 
weakly collisional at several arcseconds, the conductiv- 
ity is well approximated by a collisionless formula with 
k cx r. Another assumption of the kind deals with the 
electron velocity. As electrons can get only mildly rela- 
tivistic, we take conductivity to be proportional to square 
root of electron temperature n cx y/T^, instead of p ropor 
tionality to relativistic electron velocity K cx v c (John 
son & Quataert 2007). When the gradient of electron 



temperature gets too large, the electrons transport heat 
via a constant saturated flux, instead o f the flux propor- 
tiona l to the gradient of temperature (Cowie & McKee 
19761. We check a posteriori that the flow is in an un- 



saturated heat flux regime. Finally, we have for the heat 
flux Q = —nkBdT e /dr 



n = 0.1y/kBT e /m e rn, 



(1) 



where n = n e is the electron density (Cowie & McKee 



4 



Shchcrbakov R., Baganoff, F. 



19761. 



4.2. System of Equations 

Gravitational energy of gas in the potential of an ac- 
cretor is the ultimate inflow driver. It gets transformed 
directly in several types: kinetic energy of bulk toroidal 
and radial motion, energy of turbulent magnetic and ve- 
locity fields, thermal energy. Turbulent energy can also 
originate from the toroidal shearing flow in a disk. Tur- 
bulence dissipates into thermal motions of ions and elec- 
trons on the dynamical timescale, whereas ions and elec- 
trons exchange energy by slow Coulomb collisions. The 
faster collective modes of ion-electron energy exchange 
may exist, tho ugh they may not lead to equilibration of 
temperatures (jShkarofs ky et aL||1966 ). We do not sepa- 
rate the turbulent term or write an equation on it for the 
purpose of current work, as its direct dynamical influence 
is smaller than the influence of additional thermal en- 
ergy produc ed via dissipation o f turbulence and entropy 
productio n ( Sh cherbakovj |2008[ ) . Following Johnson &; 
Quataert (2007), we introduce the fractions f p and / e 
ot changes ot gravitational energy, which go directly into 
thermal energy of ions and el ectrons, but relate the m via 



a direct heating mechanism (Sharma et al. 2007a I . For 



the purpose of numerical stability we enhance Coulomb 
collisions by a factor of 1000, which effectively makes ion 
and electron temperatures equal at large distances from 
the BH, but does not influence T e near the BH. Let us 
convert the qualitative ideas into equations. 

The composition of plasma determines the exact bal- 
ance of the black hole gravitational pull and support- 
ing gas pressure. Let us define the source function q, 
so that the ejected mass of stellar winds per second is 
Mw = J 47r r 2 q dr. We denote the electron density by 
n = n e and write the continuity equation as 



dn 1 d(nv r r 2 ) 
dt r 2 dr 



where 



Ma 



1.14 



(2) 



(3) 



is the average atomic mass per one electro n for assumed 
solar abundance of fully ionized elements (Najarro et al. 
2004 1 . The ratio of number densities of atomic nuclei to 
electrons is 

d = ra non _ e i/n « 0.93. (4) 

We write separate energy equations for electrons (e) 
and all ions (p) in terms of 



lk B T e 



and 



lk B T p 



-"Sp 



(5) 



assuming all ions have the same temperature. We set 
the speed of light equal unity c = 1 and normalize to 
it all velocities. The ideal gas law gives normalized gas 
pressure 



Pgas = Pp + Pe = n(c 2 se + d ■ C 2 sp ) 

to be substituted into the Euler equation 
Dv r 1 <9p gas r g q(r) 



Dt 



n/j, av dr 2(r - r g ) 2 n/i av 



(6) 



0, (7) 



where D/Dt — d/dt + v r d/dr. The last term corresponds 
to zero bulk radial velocity of emitted stellar winds. 

The electron internal energy density can be approxi- 
mated as 



3^3 (e- 



■Kil 



3 0.7 + 2c 2 se m p /m e 
2 0.7 + c 2 e m p /m e 



1 



(8) 



This take s into account the differ ential heat capacity of 
particles (Shkarofsky et al. 1966). The ion internal en- 
ergy per particle is u p = 3/2m p Cg p . 

The energy exchange rate by Coulomb collisions is 
( jShkarofsky etaLp9l56 ) 

„2 



F pe = 4.3 • 10 



-19 



(c 



(9) 



The non-relativistic formula is used everywhere, as F pe 
rate is only significant in the region of non-relativistic 
electrons. The energy equation for electrons is then 

'3 0.7 + 2c 2 e TO p /m e 
2 



D 



Dt 



-f e n 



3 m p /m e 
v r q(l 



Dn 



Dt 



CK 



pe 



d) 



2r 2 



2/i Q 



10) 



+ — d r (r 2 nd r c 2 se ), 



where C ~ 1000 is the enhancement of collisions and 
conductivity is given by equation ([I]) . The left-hand side 
of the equation ( 10 ) represents the compressive heating 
in the adi abatic flow. The Paczynsk i-Wiita gravitational 
potential ( |Paczynski fc Wiita 19801 is implemented for 
gravitational force, but not in the entropy production 
term. This reflects the fact that the dissipation of tur- 
bulence ceases near the BH as having slower timescale 
compared to the inflow time. The energy equation for 
ions reads 



D 

''Dt 



3 c 2 
2 sp 



-f P n 



r g v r q(l + d) 



2r 2 



+ 



2[i a 



2 °n 
sp Dt 



= -CK 



pe 



2 Csp 



(11) 



The energy injection rate into ions is chosen to be the 
same per electron as the energy injection rate into elec- 
trons to facilitate the equality of ion and electron temper- 
atures. Let us write a condition on f p and f e to decrease 
the number of free parameters. We assume the ratio of 
heatin g fractions to be given by the direct heating mech- 
anism (Sharma et al. 2007a) as 




(12) 



despite this calculation is non-relativistic and a large 
fraction of energy dissipates at the small scales instead 
of direct large-scale heating. 

5. SOLUTIONS AND DISCUSSIONS 

We solve the derived system of equations from the 
outer boundary of the feeding region at 14" = 1.3 • 10 6 r g 
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to the inner boundary at about 1.3r g , thus covering 6 or- 
ders of magnitude in radius. Such a huge dynamic range 
requires the special solution technique, the solution of 



a time-dependent system of equations (Quataert 20041 
not being an option. We employ the snooting method 
and find the smooth transonic solution through the inner 
sonic point at ~ 3r g . In the presence of conduction the 
point, where sound speed equals inflow velocity is not 
special anymore, and instead the point, where isother- 
mal speed equals the inflow velocity, plays the role of 
transonic surface (Johnson & Quataert 2007). The sys- 
tem of equations is reduced to one temperature in the 
outfl ow b y s etting T e = T p and adding the equations 
(10) and (11). The inner boundary is set at a point n n , 
where dT e /dr — in a non-conductive solution. Then 
for any non-zero conductivity the zero heat flux condi- 
tion dT e /dr — is enforced at r- m . The outer boundary 
condition at r out is uncertain. It is natural to think the 



outflow would be transonic (Lamers & Cassinelli 1999), 
however, significant outer pressure may hold the gas m 
the subsonic regime near r ou t- The position of zero veloc- 
ity stagnation point r st determines the accretion rate M . 
Instead of setting the pressure at the outer boundary we 
regulate that pressure by setting temperature T st at the 
stagnation point. Thus, we have 4 independent variables 
in the fit: accretion rate M, temperature at stagnation 
point T st , the ion heating rate f p and the normalization 
N of the point source contribution. They are all found it- 
eratively to minimize x 2 - We also iteratively find the po- 
sitions of sonic point and stagnation point. The positions 
of inner boundary and outer boundary are unchanged 
while solving the 4-point boundary value problem. 

The observed surface brightness radial profile is the 
data we fit. We generate a surface brightness profile cor- 
responding to the dynamical model by performing the 
optically thin ray tracing of X-Rays at a set of photon 
energies and projected distances from the BH. We em 



loy t he up-to-date bremsstrahlung emissivities ( Gould 
( 1980) and errata) and account for the emission by heavy 
elements, e xcluding iron. Solar metallicity interstellar 



absorption (Morrison & McCammon 1983) is assumed 
with hydrogen column TVh = lO^cm -1,2 . The fluxes are 
convolved with the response of Chandra to find counts, 
then blurred with the energy-independent PSF (see Fig- 
ure^ and integrated over the radial extent of each bin. 
The model with M = 6 • lO^Mgyear" 1 , f p = 0.46, 

T st = 3.2 • 10 7 K and 550counts pixel -2 produced at 
r = by a point source gives an excellent fit with the 
minimum reduced \ 2 — 1.45 and weighed xtd = 0-68 
with 1/r weights. The stagnation point is at r st = 1.01". 
The correspondent unabsorbed point source luminosity 
L = 4 ■ 10 32 crg s -1 is estimated for monoenergetic pho- 
tons at 4 ke V and agrees with the estim ates of SSC lu- 



minosity in Moscibrodzka et al. (2009). Energy 4 keV 
is chosen as the energy Chandra is most sensitive to for 
assumed TVjj. The minimum reduced x 2 = 15 is achieved 
for the model without the point source. The models with 
the outer sonic point instead of finite bounding pressure 
underpredict the X-ray surface brightness at several arc- 
seconds, assuming fixed = 10 23 cm -2 . The reliable fit- 
ting for iV H is possible only with the use of spectral data 
and is left for future research. The assumption T p = T e 
represents the additional point of concern. Temperature 



equili brium might not hold at the stagnation point at 
1" ( Quataert]|2004 ) , however the thermahzation rate ex- 
ceeds the outflow rate at 5" in our subsonic dense out- 
flow, thus Tp — T e holds there. The reliable modeling of 
non-equilibrium flows requires the modeling of the whole 
spatial structure of the stellar winds and is left for the 
future research as well. 



Outflow velocity 300km/s 




-J. 



Baganoff et al.(2003): 
n e =130cm- 3 ? r ( ,=2keV at 1.5' 




Fig. 4. — Radial profiles of electron density n = n e in cm 
(upper panel) and electron temperature T e in keV (lower panel) in 
the feeding region. 



The profiles of electron density n e and temperature T e 
within several arcseconds from the BH are shown on Fig- 
ure |4land_£OTnmrej^ll wi earlier estimates 
( |Baganoff et al.||2003| |Quataert|[2004| . The difference 
is that our best fit is a subsonic flow supported by the 
outer medium with the density bounce at 5". Though the 



achieved outflow velocity v out 
dependent of radius f or r > 2" 



300km s is almost in- 



The line cooling ( Suther 



land & Dopita 1993 ) reduces the heat contents only by 



several percent for gas reaching 5", bremsstrahlung cool- 
ing being less important. 

The profiles of dimensionlcss electron tempera- 
ture fcsT e /(m e c 2 ) and ratio T p /T e within several 
Schwarzschild radii from the BH are shown on Figure [5j 
The electron temperature T e = 4 • 10 10 K and density 
n e = 2- 10 6 cm -3 are found close to the BH. This dynami- 
cal model gives an excellent fit to the optically thick lu mi- 
nosity L = 1.73 Jy at 86 GHz (|Krichbaum et al.|2006[) for 
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Inner density up to 
H e <210 6 cm -3 



Ratio T p jT e is large because of 



larger 



T p than in simulations 
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Fig. 5. — Radial profiles of dimensionless electron temperature 
normalized to electron mass kgT e / (m e c 2 ) (upper panel) and ratio 
of ion to electron temperatures T p /T e (lower panel) close to the 
BH. The inner sonic point is at 3r g . 



assumed equipartition of thermal energy with the mag- 
netic field. The model overpredicts by a factor of 20 the 
observed Faraday rotati on measure RM ~ 50cm -2 at 



230 GHz ( Marrone|2007 ), but this may well be a geomet- 
ric factor. The accretion rate, temperature and density 
near the BH are in good agreement with more compli- 
cated models specifically focusing on sub-mm emission 



(Sharma et al. 



2008) |Moscibrodzka et al||2009| ). We no 



tice that the ratio of ion and electron te mperatures T v /T e 
is s ignific antly larger than predicted by |Moscibrodzka~et| 
al. ( 2009 1 , but probably because of the significantly lower 
in their numerical simulations of the limited domain. 
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